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Abstract 

In this paper, we study the cracking of compact object PSR J1614- 
2230 in quadratic regime with electromagnetic field. For this purpose, 
we develop a general formalism to determine the cracking of charged 
compact objects. We apply local density perturbations to hydrostatic 
equilibrium equation as well as physical variables involve in the model. 

We plot the force distribution function against radius of the star with 
different parametric values of model both with and without charge. 

It is found that PSR J1614-2230 remains stable (no cracking) corre¬ 
sponding to different values of parameters when charge is zero, while 
it exhibit cracking (unstable) when charge is introduced. We conclude 
that stability region increases as amount of charge increases. 
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1 Introduction 


Self-gravitating compact objects (CO) like neutron stars, white dwarfs, mil¬ 
lisecond pulsars, and so forth, belongs to a distinguish class of those celestial 
bodies whose study become very signihcant in novel astrophysical research. 
It is evident that when a star or system of stars burns out all its nuclear 
fuel, its remnants can have one of three possibilities: white dwarfs, neutron 
stars and black holes. The stability of stellar remnants plays a key role in 
general relativity (GR) as well as modihed relativistic theories |T]. The oc¬ 
currence of gravitational collapse may be as a result of cooling of gaseous 
material, change in anisotropy, fluctuation of gravitational waves and varia¬ 
tion of electromagnetic held of CO [2] . Therefore, such phenomena stimulate 
our interest to study the stability regions of these self-gravitating CO. 

Astronomical objects are not physically viable, if they are unstable to¬ 
wards perturbations. Therefore, it is important to check the stability of these 
objects. In this context, Bondi [3] initially developed hydrostatic equilibrium 
equation to examine the stability of self-gravitating spheres. Chandrasekhar 
[1] calculated the principle value, i.e., | to determine the dynamical instabil¬ 
ity of sphere hlled with perfect huid in GR. Herrera |5] presented the tech¬ 
nique of cracking to discuss gravitational collapse of self-gravitating spherical 
CO. This technique interprets the behavior of inner fluid distribution of CO 
just after equilibrium state is disturbed. Cracking takes place in CO when 
radial forces changes its sign from positive to negative and vice versa [^. 
Several authors [Z]-[lQ] studied non-local effects of cracking through radial 
sound speed velocities and Raychaudhuri equation for spherically symmetric 
CO. Gonzalez [HI [12] presented the idea of local density perturbation (DP) 
to discussed the idea of cracking for relativistic spheres. 

To study the effect of charge on the physical properties of stars is an 
important subject in GR. In this scenario, Bonnor nanii explored the effect 
of charge on spherically symmetric GO and found that electric repulsion can 
halt the gravitational collapse. Bondi [15] used local Minkowski coordinates 
to described the contraction of radiating isotropic spherical symmetry. The 
main hindrance in astrophysics and GR is to develop stable mathematical 
models which describes the characteristic of charged spherical GO. Bekenstein 
[T6] presented the idea of gravitational collapse in charged GO. Ray et al. ini 
found the maximum amount of charge, (i.e., approximately 10^° coulomb), 
needed for GO to be in equilibrium conhguration. Some authors [HI [IH] 
studied the impact of charge on gravitational collapse of celestial objects 
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and analyzed the tendency of self-gravitating systems to prodnce charged 
black holes or naked singnlarities. Sharif and Azam [201 El] stndied the 
stability of spherical and cylindrical symmetric objects nnder the inflnence 
of electromagnetic held. 

Demorest et ah [22] nsed the Green Bank Telescope at the National Radio 
Astronomy Observatory to analyze the system of stars by means of Shapiro 
delay (SD) and presented the observed valnes of different physical parame¬ 
ters for PSR J1614-2230. These physical parameters like ecliptic longitude, 
ecliptic latitude, parallax pulsar spin, pulsar spin period, orbital period, com¬ 
panion mass, radius, and so forth are recorded with very high precision by 
SD for PSR J1614-2230. The availability of very accurate parametric values 
made PSR J1614-2230 extremely important for modern research in GR. Neu¬ 
tron stars are made of the most dense material exist in this universe. Tauris 
et ah [23] developed mathematical model of PSR J1614-2230 and provided 
the possible variation of masses to show that PSR J1614-2230 was born 
more massive as compared to any discovered neutron star. Lin et al. [23] 
used stellar evaluation code “MESA” to describe the relationship between 
PSR J1614-2230 and its stellar companion. This discovery of high massive 
neutron star has extensive consequences on the equation of state (EoS) of 
matter with high densities. The relationship between physical parameters 
become more complicated as linear EoS is replaced by nonlinear EoS. In this 
work, we apply the concept of cracking to self-gravitating GO in the pres¬ 
ence of electromagnetic held in quadratic regime. Here, we take local density 
perturbation (DP) which is different from constant DP presented by Herrera 
[5]. We applied this technique to the model of charged compact objects with 
quadratic EoS presented by Takisa et al. [25] and determine the cracking of 
newly discovered PSR J1614-2230 with electromagnetic held. Recently, we 
have investigated the cracking of some compact objects with and without 
electromagnetic held in linear regime [26] . 

This paper is arranged as follows. Section 2 deals with Einsten-Maxwell 
held and TolmanOppenheimerVolkoh (TOV) equations corresponding to an¬ 
isotropic huid. We present the general formalism to determine the cracking 
of charged GO with local DP in the quadratic regime in section 3 . Section 4 
investigate stable and unstable regions of compact star PSR J1614-2230. In 
the last section, we conclude our results. 
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2 Einstein-Maxwell Field and TolmanOppen- 
heimerVolkoff Equations 


We consider the line element for a static spherically symmetric space time in 
cnrvatnre coordinates given by 

+ r‘^{d9‘^ + sin^ 6d(j)‘^), (1) 

where 0<d<7r, O<0<27r and z/ = z^(r), A = A(r) are gravitational 
potentials. The Maxwell’s eqnations are dehned as 


F ab\c FjyQ-a -fca;6 0, 

(2) 

4nr, Eab = FaA^ - \gabFadF^^, 

(3) 


where is the electromagnetic held tensor, J is the fonr cnrrent density 
and Eab is the electromagnetic energy-momentnm tensor [27]. The skew- 
symmetric electromagnetic held tensor can be decomposed as 
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where E = [E^, Ey, E^) is the electric held and B = {B^-, By, Bz) is the 
magnetic held. The electromagnetic held tensor and fonr cnrrent density can 
be dehned as 

Fab = A,a - A,b, A = (yA, (5) 

where A and a are the fonr potential and proper charge density and = 
is fonr vector velocity of the hnid. The fonr potential is dehned as 


/l„ = Wr), 0,0,0). 


( 6 ) 


Using this in above eqnation, it yields 

^01 = 

which can also be written as 

FOi = 


(7) 
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( 8 ) 




where, we have used E{r) = e The total energy-momentum ten¬ 

sor corresponding to charged anisotropic fluid sphere is deflned by [Tf] 

j^2 j^2 j^2 

Tab = diag{—p - 

The terms E, p, Pr and Pt are electromagnetic held, energy density, radial 
and tangential pressure respectively. 

The synergies of electromagnetic held and matter are governed by system 
of held equations. These synergies of spherically symmetric metric corre¬ 
sponds to Einstein-Maxwell held equations given by 

Gab = iGTab = n{Mab + Eab) (10) 

where Mab is the energy momentum tensor for the fluid inside the star and 
Eab = FacF^ — jdatFcdF'^'^ is electromagnetic held tensor. The non-zero 
components of Einstein-Maxwell held equations corresponding to Eqs. ([T]) and 
do]) are given as follows 


1 + e“2^(2rA' - 1) 

2 2^ 

= diTTrypP — , 

(11) 

1 - e“2^(2rz/' + 1) 

2 2^^ 

= Pr + , 

(12) 

+ u'X'r^ - 

2 2^^ 

= -Svrr^Pi-r^ —, 

( 13 ) 

r‘^a 

= e-\r^E)\ 

( 14 ) 


where “/” denotes the differentiation with respect to r. 

It is clear that the choice of EoS of fluid inside the star plays a key role 
for its physical significance. Thus, a star is physically acceptable, if it satisfy 
the barotropic EoS Pr = Prip)- In this work, we have used the quadratic 
EoS to explore the stability of PSR J1614-2230. The quadratic EoS is given 
by m 

Pr = If? + ap- j3, (15) 

where 7 , a and {3 are constants and are constrained by (p < 1^) and 
(3 = ape, where pe = 0.5 x IQ^^g/crn? gives the density at the boundary 
surface of sphere. It is interesting to note that this equation reduce to linear 
EoS, when 7 = 0 [25] . 
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Solving Eqs. (fTT]) - (fT^ simultaneously, we obtain hydrostatic equilibrium 
equation (TOV) for anisotropic charged fluid 


dPr 2{Pt - Pr 


E 


(p + P,)„' + —(r^EY, 


( 16 ) 


dr r \r • ' / • 

which shows that gradient of pressure is effected by charge and anisotropy of 
fluid. Using the relation = 1 — 2M/r + Q'^/r'^ in the above equation 

j, it yields 


^ dPr 2(Pt-Pr) , + Ar^E^ -STTT^Pr (r‘^E)'E , , 

f! = —i; + -Lj —+ (p + p,) + -Pp- = 0.(17) 


dr ' r ' \r • I / _ 2^ _j_ ^ 2 ^ 2 ^ 

where the mass function with Q = r‘^E is defined as 


dvrr^ 


M = 4:71 {p{x) + —)x^dx. 

. n OTT 


(18) 


3 Effect of Local Density Perturbation 


In this section, we perturb the equilibrium conhguration of charged CO 
through local DP {6p). Eq. (lT7)l depicts that cracking take place in interior 
of spherical CO when equilibrium state is interrupted due to change in sign 
of perturb force i.e., <50 < 0 —)■ hO > 0 and vice-versa. We apply the local 
DP to Eq. dTT)) and all the physical variables like mass, radial and tangential 
pressure, electromagnetic field and their derivatives involve in Eq. fll7p . given 
by 


Pr{p + 6p) = P,{p) + ^6p, (19) 

dPj. I r \ dPp /X d f dPp\ dPp d p \ /onN 

_(p + ip) = (20) 

dP 

Pt{p + 5p) = Pt{p) + -j^5p, (21) 

M{p + 5p) = M(p) + ^5p, (22) 

dp 

E{p + dp) = E{p) + (23) 

E\p + 5p) = E\p) + ^5p. (24) 
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The radial sound speed and tangential sound speed are define as 
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The perturb form of Eq. (fT7|) is given by 

= f]o(p, Pr, p;, Pt, M, E, E') + 

where 
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(28) 


This is the fundamental equation used to determine the effects of local DP on 
the cracking of charged anisotropic fluid. We will plot the force distribution 
function against radius “r” of the star for different values of the parameters 
involve in the model. Using Eq. fITT)) . the derivatives involve in the above 
equation are given as follows 
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Figure 1: Plots shows that there is no cracking, i.e., the PSR J1614-2230 
remains stable for different values of the parameters involved in the model 
given in Table 1, when = 0 in quadratic regime. 


4 Cracking of PSR J1614-2230 


Here, we apply the formalism developed in the above section to investigate 
the cracking of charged objects for the model given by Takisa et al. |2S|. 
This model is consistent with the physical features of observed objects and 
its connection can be made with PSR J1614-2230 for particular values of 
parameters given in [25] . The analysis of Takisa was seems to be consistent 
with observational objects such as Vela X-1, Cen X-3, SMC X-1, PSR J1903- 
327 and PSR J1614-2230. But our focus in this analysis is the particular 
object PSR J1614-2230 because its mass and radius has been measured with 
great accuracy. The model is dehned by following equations 
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8(ar2 + l) + 
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Figure 2: Cracking of PSR J1614-2230 with 7 = 0.0, a = 0.99 and si = 
0,7.5,14.5. 
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The radial and tangential sonnd speed velocities can be obtained from Eqs. fl37)l 
and fl38|) as 

= a + 27 p, (43) 

and 
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Figure 3: Cracking of PSR J1614-2230 with 7 = 0.140, a = 0.33 and Si 
0,7.5,14.5. 



Figure 4: Cracking of PSR J1614-2230 with 7 = 0.158, a = 0.24 and si 
0,7.5,14.5. 
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Figure 5: Cracking of PSR J1614-2230 with 7 = 0.163, a = 0.21 and Si 
0,7.5,14.5. 



Figure 6 : Cracking of PSR J1614-2230 with 7 = 0.177, a = 0.15 and Si 
0,7.5,14.5. 
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Figure 7: Stability regions for 7 = 0.196, a = 0.06 and Si = 0, 7.5,14.5 
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Figure 8 : Stability regions for 7 = 0.200, a = 0.04 and si = 0, 7.5,14.5 
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Table 1: Stability of neutral PSR J1614-2230 when 7 and hi are variable 


7 

ai 

hi 

a 

r(km) 

Rc(km) 

0.100 

53.34 

6.90 

0.33 

11.07 

Stable 

0.126 

53.34 

8.74 

0.33 

10.85 

Stable 

0.132 

53.34 

10.74 

0.33 

10.60 

Stable 

0.140 

53.34 

13.33 

0.33 

10.30 

Stable 

0.148 

53.34 

15.61 

0.33 

9.99 

Stable 

0.154 

53.34 

16.87 

0.33 

9.82 

Stable 

0.163 

53.34 

19.04 

0.33 

9.51 

Stable 

0.177 

53.34 

21.72 

0.33 

9.13 

Stable 

0.189 

53.34 

23.64 

0.33 

8.83 

Stable 

0.196 

53.34 

24.73 

0.33 

8.65 

Stable 

0.200 

53.34 

28.42 

0.33 

8.04 

Stable 


Table 2: Stability of PSR J1614-2230 when si = 0 and a, 7 are variable 


7 

Oi 

hi 

a 

r(km) 

Rc(km) 

0.0 

53.34 

13.33 

0.99 

10.30 

7.7 

0.140 

53.34 

13.33 

0.33 

10.30 

Stable 

0.158 

53.34 

13.33 

0.24 

10.50 

Stable 

0.163 

53.34 

13.33 

0.21 

10.70 

Stable 

0.177 

53.34 

13.33 

0.15 

10.90 

Stable 

0.196 

53.34 

13.33 

0.06 

11.06 

Stable 

0.200 

53.34 

13.33 

0.04 

11.09 

Stable 


The constants a, b and s have dimension of length (L and chosen in 
such a way that the given system satisfy the following conditions 

• energy density must remains positive before and after equilibrium state. 
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Table 3: Stability of PSR J1614-2230 when si = 7.5 and a, 7 are variable 


7 

ai 

bi 

a 

r(km) 

Pc (km) 

0.0 

53.34 

13.33 

0.99 

9.67 

8.6 

0.140 

53.34 

13.33 

0.33 

9.67 

8.6 

0.158 

53.34 

13.33 

0.24 

10.07 

9.3 

0.163 

53.34 

13.33 

0.21 

10.37 

9.3 

0.177 

53.34 

13.33 

0.15 

10.56 

9.5 

0.196 

53.34 

13.33 

0.06 

10.65 

9.6 

0.200 

53.34 

13.33 

0.04 

10.64 

9.7 


Table 4: Stability of PSR J1614-2230 when si = 14.5 and a, 7 are variable 


7 

di 

bi 

a 

r(km) 

Pc (km) 

0.0 

53.34 

13.33 

0.99 

9.21 

8.3 

0.140 

53.34 

13.33 

0.33 

9.21 

8.3 

0.158 

53.34 

13.33 

0.24 

10.05 

9.1 

0.163 

53.34 

13.33 

0.21 

10.10 

9.1 

0.177 

53.34 

13.33 

0.15 

10.15 

9.2 

0.196 

53.34 

13.33 

0.06 

10.18 

9.2 

0.200 

53.34 

13.33 

0.04 

10.19 

9.2 


• radial pressure should be vanishes at the boundary of star. 

• At the center of star i.e. r = 0, we have = A = 0. 

• is constant in the quadratic regime. 

• Across boundary of star, when r = e, we have 
g-2A _ ^ _ 2M/e + Q^je^ 

^ ^ _ 2Mle + 

By considering above conditions, we have 


a = 0.33, 




where 3fJ=43.245 km and the values given above 
vational values given by Takisa et ah 


Si 

are compatible with obser- 
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For the sake of regions (stable and nnstable) of PSR J1614-2230, we have 
plotted force distribntion fnnction against radius of the star for different 
values of the parameters involve in the model shown in Figures 1 - 8 . We 
summarizes these results as follows: 

• Figure 1 depicts that all curves do not change its sign with different 
values of 7 and 61 corresponding to Table 1 . Hence, we hnd that PSR 
J1614-2230 is stable in the absence of charge in quadratic regime and 
it is unstable in linear regime which is analogous to the results found 
in [SB]. From Table 2 , it is clear that any variation in coefficients of 
quadratic EoS does not effect stability even radius of PSR J1614-2230 
changes approximately to 4%. 

• In Figure 2 , there are three curves corresponding to model parameters 
7 = 0.0, a = 0.99 and charge si = 0, 7.5, 14.5. It is noted that all 
three curves (red, blue and green) change its sign for charge parameter 
(si = 0, 7.5, 14.5), respectively in the linear regime. This shows 
that PSR J1614-2230 is unstable in linear regime, where the symbols 
“ o ”, “o”, “ * ” represents the cracking points (where curve changes 
its sign from negative to positive) corresponding to Si = 0, 7.5, 14.5, 
respectively. The cracking values Rc = 7.7, 8 . 6 , 8.3 corresponds to 
Si = 0, 7.5, 14.5 (red, blue and green) are given in (Table 2 , 3 , 4 ). In 
this case, our results are consisted with [26] in linear regime. 

• Figures 3-8 represents the cracking of PSR J1614-2230 star for hxed 
values of parameters 7 = 0.140,0.158,0.163,0.177,0.196,0.200, a = 
0.33, 0.24, 0.21, 0.15, 0.06, 0.04 and charge Si = 0.0, 7.5,14.5 in quadratic 
regime. We see that cracking take place for charge parameter si = 7.5 
and 14.5, which are indicated by the cracking points “o” and “ * ” 
corresponding to blue and green curves, respectively. These cracking 
points {Rc) are given in Table 3 - 4 . However, in each case, the star re¬ 
main stable, i.e., no cracking take place for si = 0 in quadratic regime. 
Hence, PSR J1614-2230 exhibit cracking both in linear and quadratic 
regime when charge is present. From these illustrations, we conclude 
that as charge increases cracking points are slightly shifted towards 
center, which indicates that binding forces of CO become stronger and 
more mass is directed towards origin. 
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5 Conclusions and observations 


We have applied the technique of cracking presented by Herrera [5] to charged 
anisotropic self-gravitating CO. The impact of local DP on the stability of 
inner fluid of star in the presence of charge is considered in the scenario of GR. 
It has been observed that cracking of CO takes place when the system leave 
its equilibrium state. The numerical value of Rc (cracking point) provides 
the stable/unstable region in the quadratic regime. 

We have used the model of Takisa et ah [25] to investigate the cracking of 
PSR J1614-2230 with and without charge. Figure 1 represents the stability 
of celestial object PSR J1614-2230 in the absence of charge for different 
values of parameters 7 and hi given in Table 1 . It is shown that PSR J1614- 
2230 remains stable, when quadratic EoS is considered in neutral case but it 
exhibited cracking with the inclusion of charge. Figure 2 has been plotted 
for 7 = 0.0 and different values of a and charge Si = 0, 7.5, 14.5. It is shown 
that PSR J1614-2230 exhibit cracking in each case given by Rc = 7.7, 8 . 6 , 8.3. 
For Si = 0, PSR J1614-2230 shows cracking which is consistent with our 
recent published work [26]. It is worth mentioned here that our results are 
analogous to [26], when 7 = 0 (Linear Regime) in the presence of charge. 

In hgures 3-8, we have given the comparison of stability region with dif¬ 
ferent values of a, 7 and charge parameter Si. In these hgures stability 
regions are plotted for values of charge parameter si = 0, 7.5,14.5 with Rc 
represented by “ o ”, “o” and “ * ” for si = 0.0, si = 7.5 and si = 14.5, re¬ 
spectively. We observe that the value of Rc decreases as electromagnetic held 
increases, which are given in Tables 2-4 for diherent values of parameters. 
Figures 3-8 shows that cracking take place in each case for diherent values 
of the parameters corresponding to Si = 7.5 and Si = 14.5 in the quadratic 
regime, while remains stable in the absence of charge (si = 0 . 0 ). 

It is noted that the local DP scheme does not ahect the stability of CO 
(remains stable) in neutral case, while change its stability (potentially un¬ 
stable) drastically with the inclusion of charge in quadratic regime. Thus, 
the local DP scheme under nonlinear EoS considerably ehect the stability 
regions of CO. When physical parameters like mass, electromagnetic held 
and density of anisotropic charged self-gravitating CO are locally perturbed, 
they drastically ahect the sensitivity of radial forces which may leads towards 
the gravitational collapse. Therefore, the stability region of PSR J1614-2230 
increases as value of electromagnetic held increases. Hence, we conclude that 
the binding forces of neutron star PSR J1614-2230 become stronger as we 
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move towards center of star and it becomes more dense as charge increases. 

It is important to mentioned here that the idea of cracking was presented 
by Herrera |5] to understand the behavior of inner fluid distribution just after 
departure from equilibrium state may be responsible for cracking (overturn¬ 
ing) of anisotropic sphere |8]. In his study, the global DP affects physical 
quantities like mass, tangential and radial pressure but do not effect pressure 
gradient. In this work global DP technique is modihed by local density per¬ 
turbations to study cracking in the presence of electromagnetic held. Finally, 
we conclude that the the given object exhibits cracking in the presence of 
electromagnetic held. 
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